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■ We employ a helicity-basis kinetic equation approach to investigate the anomalous Hall effect in 

(N ■ two-dimensional narrow-band semiconductors considering both Rashba and extrinsic spin-orbit (SO) 

couplings, as well as a SO coupling directly induced by an external driving electric field. Taking 
^ . account of long-range electron- impurity scattering up to the second Born approximation, we find that 

' the various components of the anomalous Hall current fit into two classes: (a) side-jump and (b) skew 

scattering anomalous Hall currents. The side-jump anomalous Hall current involves contributions 
not only from the extrinsic SO coupling but also from the SO coupling due to the driving electric 
field. It also contains a component which arises from the Rashba SO coupling and relates to the 
off-diagonal elements of the helicity-basis distribution function. The skew scattering anomalous Hall 
, effect arises from the anisotropy of the diagonal elements of the distribution function and it is a 

^ ■ result of both the Rashba and extrinsic SO interactions. Further, we perform a numerical calculation 

to study the anomalous Hall effect in a typical InSb/AlInSb quantum well. The dependencies of the 
side-jump and skew scattering anomalous Hall conductivities on magnetization and on the Rashba 
SO coupling constant are examined. 
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I. INTRODUCTION 



Q , A nonvanishing magnetization in spin-split systems may lead to an extraordinary Hall current.^ This so-called 
Q • anomalous Hall effect (AHE) was first observed more than a century ago,^ but its complete understanding still 
I— remains a challenge today. Up to now, it has been made clear that the AHE arises essentially from spin-orbit (SO) 
^ • interactions, which are the results of rapid movements of carriers in various electric fields, such as nuclear fields {e.g. 
the Dresselhaus SO coupling), electric fields associated with strains or gate biases [e.g. the Rashba SO coupling), or 
. fields induced by electron-impurity scattering (extrinsic SO coupling), etc^ 
,—1 ■ In 1954, Karplus and Luttinger proposed, for the first time, a mechanism of the anomalous Hall effecti^ This 
, mechanism is associated with the spin-orbit interaction due to nuclear fields and yields to an anomalous Hall current 
0^ (AHC) independent of any electron-impurity scattering. Recently, it has been reformulated by Jungwirth et al. within 
, a framework of momentum-space Berry phase,® and was used to explain the AHE in various ferromagnetic systems, 
such as dilute magnetic semiconductors (Ga,Mn)As,®- ferromagnetic Fe^ SrRuO series ^iiSiii spinel CuCrSe etci^ 
The AHE may also stem from a spin-orbit coupling induced by electron-impurity scattering, i.e. the extrinsic 
"j^ SO coupling."^ It was found that there are two mechanisms responsible for this extrinsic AHE: a side-jump process 
' proposed by Berger— and a skew scattering given by Smitii^ The side-jump AHE arises from a sidewise shift of the 
'J^ center of the electron wave packet and relates to an anomalous term in the current operator caused by the extrinsic SO 
coupling. ^^'^5 The skew scattering AHE corresponds to an anisotropic enhancement of the wave packet due to electron- 
^ , impurity scattering and can be accounted for by considering the scattering in the second Born approximation^ It 
O ' was also clear that the side-jump anomalous Hall conductivity is independent of impurity density nt, while the skew- 
, scattering one is proportional to (n^)^^. Recently, considering a short-range electron-impurity scattering, Crepieux, 
^ I et al. presented a unified derivation of both the side-jump and skew scattering mechanisms within the framework of 
' a formal Dirac equation for the electrons. The weak- localization corrections to these anomalous Hall currents have 
, also been investigatediiSiiSiSfliSi 
d ' The Rashba SO coupling in a two-dimensional (2D) electron system with magnetization can also give rise to a 
nonvanishing contribution to the Hall conductivity^^'^^'^'* It was found that the anomalous Hall current due to the 
Rashba SO interaction consists of two terms: one of which is associated with all electron states below the Fermi 
surface and is independent of any electron-impurity scattering; and another one relating only to electrons near the 
Fermi surface, which is disorder-related but independent of impurity density^2*24 Liu and Lei have also clarified 
that, in the helicity basis, these two different contributions to anomalous Hall current in 2D Rashba semiconductors 
essentially arise from two distinct interband polarizations and they relate to two distinct parts of the off-diagonal 
elements of the distribution function^ In these studies, the Rashba SO coupling was considered nonperturbatively, 
but the extrinsic SO interaction was completely ignored. Such a treatment is valid only for 2D electron systems based 
on wide-band semiconductors. In 2D narrow-band semiconductors, such as InSb/AlInSb quantum wells (QWs), the 
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coupling constant for extrinsic SO interaction is relatively large (for example, the extrinsic SO coupling constant, A, 
is A = 5.31 nm^ for InSb, while it is equal to 0.053 nm^ for GaAs^^). Hence, to investigate the AHE in Rashba 2D 
systems based on narrow-band semiconductors, we have to consider not only the Rashba but also the extrinsic SO 
interactions. 

In this paper, we employ a kinetic equation approach to investigate the AHE in Rashba 2D narrow-band semi- 
conductors. We deal with the Rashba SO interaction in a nonperturbative way, while the extrinsic SO coupling is 
considered in the first order of the coupling constant A. We also take account of the SO interaction induced directly 
by the external driving electric field, which, to our knowledge, was mentioned only by Nozieres and Lewiner in the 
absence of Rashba SO couplingi2& In our study, to investigate the skew-scattering AHE efi'ect, we consider the electron- 
impurity scattering up to the second-Born approximation. It is found that various components of the anomalous Hall 
current can fit into two classes: the side-jump and skew scattering anomalous Hall currents. The side-jump AHC 
involves contributions from the extrinsic SO coupling and SO coupling induced by the driving electric field. It also 
contains a component which comes from the Rashba SO interaction and relates to the off-diagonal elements of the 
helicity-basis distribution function. The skew scattering AHC is associated with anisotropic diagonal components of 
the distribution function and stems from both the Rashba and extrinsic spin-orbit interactions. A numerical calcu- 
lation of the anomalous Hall current in a InSb/AlInSb quantum well indicates that both the side-jump and skew 
scattering anomalous Hall currents are of the same order of magnitude, leading to complicated dependencies of the 
total anomalous Hall conductivity on magnetization and on the Rashba spin-orbit coupling constant. It is also clear 
that in the side-jump anomalous Hall current, the contribution from SO coupling due to the driving electric field is 
dominant for small magnetization. 

This paper is organized as follows. In Sec. II, we derive the kinetic equation for the nonequilibrium distribution 
function considering long-range electron-impurity scattering up to the second Born approximation. The solution of 
this equation to first order of the extrinsic spin-orbit coupling is presented. We also discuss the various components 
of the side-jump and skew scattering anomalous Hall currents. In Sec. HI, we perform a numerical calculation to 
investigate the anomalous Hall effect in a InSb/AlInSb quantum well. Finally, we review our results in Sec. IV. The 
detailed form of the scattering term of the kinetic equation is presented in an Appendix. 



II. FORMALISM 



A. Hamiltonian and current operator 



We consider a Rashba quasi-two-dimensional narrow-band semiconductor in the x — y plane. When a homogeneous 
magnetization M = (0,0, A/) induced by a weak magnetic field B, M — g^sB {g is the effective g-factor and 
is the Bohr magneton), and a uniform in- plane dc electric field E are present, the Hamiltonian of an electron with 
momentum p = {px,Py) = ipcos(j)p,psm(/)p) can be written as 

H = Ho + i/inip + He- (1) 

Hq is the free electron Hamiltonian given by 

iJo = £p + a{<y''Py - (jypx) - Ma\ (2) 

where a is the Rashba SO coupling constant, a' {I = x, y, z) are the Pauli matrices, and £p — p^ /2m* with m* as the 
electron effective mass. By a local unitary spinor transformation, C/p, 
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" yf2\p V Je'^p^V^ -je**p^Ap + Af ) ' ^''^ 

Hamiltonian 101 can be diagonalized as Hq = UpHoUp = diag(ei(p), S2{p)) with ep(p) = p^ /2m* -t- (— l)^Ap (/x = 1, 2 
and Ap = Af ^ + a^p"^) as dispersion relations of two spin-orbit-coupled bands. 

Since the extrinsic spin-orbit coupling in narrow-band semiconductors cannot be ignored, the Hamiltonian -ffimp, 
which describes the electron-impurity interaction, should contain not only an ordinary scattering potential term but 
also a term related to the extrinsic SO coupling: 



where r and R,;, respectively, are the coordinates of the electron and impurity, V{r) is the electron-impurity scattering 
potential, and A is a spin-orbit coupling constant depending on the intrinsic semiconductor parameters, such as energy 
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gap i?Oi spin-orbit splitting Asoi and matrix element of the momentum operator between the conduction and valence 
bands P: A = [1/Eq — 1/{Eq + Aso)^]P^/3j^ He describes the application of the external electric field, and, in the 
Coulomb gauge, it can be written as 

He = -eE • r - A[cr x E] ■ p. (5) 

In Eq. ISJ, we have considered the effect of the spin-orbit coupling directly induced by the external driving dc electric 
field. 

From Hamiltonian it follows that, in spin basis, the single-particle current operator, j(p), can be written as 

i/(p)=jf(p)+ir(p)+if(p), (6) 

with I = x,y. The term j;(p) comes from the free-electron Hamiltonian Hq: j[(p) — epi/m* — aeimz'^"^ {fn = x,y 
and eimz is the totally antisymmetric tensor), while j™^{p) comes from the SO coupling term of -ffimp and takes the 
form (n ~ x,y) 

i7P(p) = iAe^ yp_ke'^'-('^-P) [£/,„„ - p„,)<j"]. (7) 

k,i 

The term jf'(p) arises from the spin-orbit coupling directly induced by the external driving electric field and it is 
given by 

if (p) = -Xehimna^'En. (8) 
Taking the statistical ensemble average, we find that the observed net current, J, consists of three components: 

Ji = Ji + j'r + jF- (9) 

Ji'^ is determined by jj'^ — Tr[jJ'^(p)/5(p)], with p{p) as the distribution function in the spin basis: p^i/(p) =< 
Ptpi'np > (V'mp ^^'^ V'a'p> respectively 

from the current operator term j™''(p) and takes the form 



''Pt-p'4'nY> > (''/'mp and 'i/'^p, respectively, are the spin-basis electron creation and annihilation operators). J™'' arises 



jr^ = ixe J2 Vp_ke'^-(''"p'{<v;+p4k> [ei™„(fc™-p,n)(7::j}. (lo) 

p,k,z,/i,i/ 

Obviously, to determine J""'', one has to analyze the function < tptp''P^^k >■ 

Without loss of generality, we study here the anomalous Hall current flow along the x-axis when the electric field 
is applied along the y-direction, i.e. E = {0,E,0). In helicity basis, the current jf can be written as 

= \e^EY,[Piiip)-P22ip)], (11) 
p 

with Ppiy(p) (fJ-jV = 1,2) defined as the elements of the helicity-basis distribution function related to the spin-basis 
distribution function by p{p) = Up p{p)Up. can be expressed as a sum of the contributions from the diagonal and 
off-diagonal elements of the helicity-basis distribution function, and J^°: 

Jx ~ Jx + "^x I (12) 

with 

= I {^-^ - pcos(/)p/5ii(p) -I- + pcos0pP22(p)| , (13) 



p 

and 



f° = |^^^cos0pRe/5i2(p) -f 2asin0plm/5i2(p)| . (14) 
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From Eq. (|14|l . we can see that the nonvanishing contribution to comes from the component of the real part of 
Puip), depending on the momentum angle through coS(/)p: Re[(/5)i2(p)] — £,{p) cos0p + and from the component 
of the imaginary part of /5i2(p), involving sin(/)p: Im[(yo)i2(p)] = C{p) sin^p + .... As a result, can be rewritten as 



p 



- COS^ V.p 
A 



cos^ (j)p + C{p) sin^ 0p 



(15) 



Note that these components, J^, J""^, and , can fit into two classes: (a) side-jump and (b) skew-scattering 
anomalous Hall currents. We know that the side-jump AHE originates from the shift of the electron wave-packet 
center towards the direction transverse to the driving electric field. Such a shift is reflected by those current operator 
components involving the antisymmetric tensor eimn- fi"^^ , and the term in j[ associated with the Rashba spin- 
orbit coupling. Correspondingly, the observed side-jump anomalous Hall current J^j is a sum of J™^, and J^°: 
J|j = 4° + J™P + JE_ From Eq. JTSJ, we see that the remaining component, J^'^, becomes nonvanishing if there 
exists a component of p^^(p) depending on the angle of momentum through coS(/)p. This implies that J^J^ results 
from an anisotropy due to electron- impurity scattering and hence it just is a component of J|^: J^^ = J^. 

B. Kinetic equation and its solution 

Obviously, in order to carry out the calculation of the anomalous Hall current, it is necessary to determine the 
electron distribution function. Under homogeneous and steady-state conditions (averaging over a uniform impurity 
distribution), the helicity-basis distribution function, p(p), obeys a kinetic equation written in the form 

eE • Vpp(p) - eE • [/5(p), U+VpUp] - *AeE • (p x n)[U+a'Up, p{p)] + i[H„, p(p)] = -/, (16) 

where n is a unit vector along z-axis and / is a scattering term determined by 

/ ^Plp,c.)G<(p,c.) + E<(p,c.)G'^(p,c.)-G'-(p,c.)E<(p,c.)--G<(p,c.)snP,^)]. (17) 

G''''''^ (p, (jj) and S'''"'<(p, cj) are, respectively, the helicity-basis nonequilibrium Green's functions and self-energies. 
Eq. H16|l is derived from the Dyson equation of the spin-basis nonequilibrium lesser Green's function by applying the 
unitary transformation Up. 

In Eq. I|16|l . the electron-impurity scattering is embedded in the self-energies, Y7'°-''^{'p,uj). This interaction in the 
helicity basis is described by a potential, Vpk, which can be written as 

Vpk = U+V{p - k)C/k + iXU+[n ■ (p X k)]^^(p - k)[/k. (18) 

In terms of the Feynman diagrams, Vpk is denoted by two different interaction vertices: ordinary and anomalous 
vertices, which, respectively, are depicted in Figs. 1(a) and 1(b). Since electron-impurity scattering will be considered 
up to the second Born approximation, it is convenient to express the self-energies by means of generalized T-matrices, 
rpj^(w), which obey the equations^!. 

T'pTH = Vpk + E^pqG'-'"(q'^)^qkH- (19) 

q 

[These equations are exhibited in terms of Feynman diagrams in Fig. 1(c).] Thus, the lesser self-energy can be written 
asSa 

E<(p,c.) = n,EfpV(a;)G<(k,c.)fk"p(u;), (20) 



k 

while the retarded and advanced self-energies take the forms 

S'-'«(p,c.) = n,f;p"(a;). (21) 

In present paper, we restrict our considerations to the linear response regime. In connection with this, all the 
functions, such as the nonequilibrium Green's functions, self-energies, and distribution function, can be expressed as 
sums of two terms: A = Aq + Ai, with A representing the Green's functions, self-energies or distribution function. Aq 



5 



and Ai , respectively, are the unperturbed part and the linear electric field part of A. In this way, the kinetic equation 
for the linear electric field part of the distribution, pi(jp), can be written as 

eE ■ Vppo(p) - eE • [po(p), U+VpUp] - ^AeE ■ (p x n)[U+a'Up, poip)] + ^[Ho, Pi(p)] - (22) 

with /(^) as the linear electric field part of the collision term /: 

I^'^^ J ^[sf(p,c.)Gg(p,c.)-G<(p,c.)Eg(p,c.) + SS(p,c.)G<(p,c.)-GS(p,c.)S<(p,c.)]. (23) 

We note that, here, the effect of Gi'°'{p,uj) on distribution function has been ignored because these linear electric 
parts of the retarded and advanced Green's functions lead to a collisional broadening effect on pi{p), which plays a 
secondary role in transport studies. 

To further simplify Eq. (|22ll . we employ a two-band generalized Kadanoff-Baym ansatz (GKBA)iS2i^ This ansatz, 
which expresses the lesser Green's function through the Wigner distribution function, has been proven sufRciently 
accurate to analyze transport and optical properties in semiconductors To first order in the dc field strength, the 
GKBA reads, 

G<(p,c.) = -GS(p,c.)pi(p)-|-pi(p)Gg(p,a;), (24) 

where the retarded and advanced Green's functions in helicity basis are diagonal matrices: Gq°(p,cj) = diag[(aj — 
ei{p) ± iS)~^,{uj — £2(p) ± iS)'-^]. Note that the helicity-basis equilibrium distribution is also diagonal, po{p) — 
diag[nF(£i(p)), nF(£2(p))] with nF{uj) as the Fermi function. 

Since our studies are concerned with electron transport within the diffusive regime, it is sufficient to study the 
lowest order of the distribution function in the impurity-density expansion. From the diagonal parts of Eq. H22|l . we 
see that the leading order of the diagonal pi{p) elements, (/5i)^^(p), is proportional to {ni)~^. These diagonal elements 
of the distribution function give rise to off-diagonal elements of the scattering term which are independent of 
the impurity density. From the fact that the left-hand side of the off-diagonal parts of Eq. H22|) involves the term 
i[Ho, pi{p)] proportional to the off-diagonal elements of the distribution function: 

it follows that the leading order of the off-diagonal elements of pi{p) should be of order (rii)", i.e. independent of 
the impurity density. Note that the contributions to /'^^■' from such off-diagonal elements of pi{p) are linear in the 
impurity density and hence can be ignored, while the contributions to I'-^-' from the diagonal elements, (/5i)p^(p), 
are independent of and become dominant. Thus, to the lowest-order of the n^-power expansion, J^^^ effectively 
involves only the diagonal elements of the distribution function. 

From Eq. 122|l we see that the driving force of the kinetic equation can be classified into two classes: diagonal 
eE • VpPo, and off-diagonal — eE • [po{p), UpVpUp] and — iAeE • (p x n)[J7p a^t/p, po{p)]- In connection with this, we 
may formally split the kinetic equation into two equations with pi{p) = pi{p) + Pi^(p) as 

eE • Vppo(p) + ^[Ho, pUp)] = (26) 

- eE • [po(p), C/+Vp[/p] - *AeE • (p x n)[U+a'Up, /}o(p)] + *[i?o, ^^(p)] = 0. (27) 

From Eq. (|27|1 it is evident that /5{^(p) has null diagonal elements. Since I^^^^ depends only on the diagonal elements 
of the distribution function, p{{p) and p{^{p) can be approximately determined independently of one another by 
Eqs. ^ and lf77|) . 

Substituting the explicit forms of GQ''(p,Ci;) into Eg. ()23|l and considering Eg. H24|l . the elements of the linear 
electric-field scattering term, can be written as 

= -2Im [{f;p)^,{s^{p))] (p()^^(p) 

-27r^ 5[e,,(p) - £^(fc)](rp\);.;.(£;.(fc))(fk"p)^^(£^(fc))(p()^^(k) 

k 

-27r5^%,,(p)-£p(fc)](fp\)^p(£p(fc))(fk"p)p^(£p(fc))(p()p^(k), (28) 
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FIG. 1: Feynman diagrams for electron-impurity scattering, (a) and (b) show the interaction vertices, which, respectively, 
correspond to the original scattering potential and the potential due to the extrinsic SO coupling, (c) is the equation for 
T-matrix. (d) and (e), respectively, are Feynman diagrams for the T- matrix in the first and second Born approximations up 
to the first order of the spin-orbit coupling constant, A. 



and 

- *(r;pW(£A(p))(p()AM(p)-*(^'p"pW(£M(p))(/5(W(p) (29) 

with = 1, 2 and ft = 3 — ^. 

Further, as in all previous studies, we consider the anomalous Hall current only to the first order of the spin- 
orbit coupling constant A. Thus, the scattering term /^^^ and hence the T-matrix may be considered only in the 
lowest- and first-order of A. On the other hand, we will evaluate the diagonal elements of /^^^ up to the second-Born 
approximation, but its off-diagonal elements only in the first-Born approximation. It is widely accepted that, in usual 
cases, the self-consistent first Born approximation is sufficiently accurate to analyze transport in diffusive regime 
(correspondingly, the scattering term may be considered only in the first-Born approximation). In our studies, ll]} is 
evaluated up to the second Born approximation because we should account for the skew scattering AHC associated 
only with the diagonal elements of the distribution function. Under these considerations, we need to analyze only 
the Feynman diagrams depicted in Figs. 1(a), 1(b), 1(d), and 1(e). Substituting the T-matrix obtained from these 
diagrams, we obtain the explicit form of the scattering term which is presented in the Appendix. 

Considering the elastic features of the electron-impurity scattering, Eq. (|26|) can be solved analytically. We know that 
/(^^ does not involve the off-diagonal elements of /5{(p). Hence, the diagonal /5{(p) elements can be determined from the 
diagonal parts of Eq. H26|) . while its off-diagonal elements are given by substituting the obtained results for {p{)^fj.{p) 
into the off-diagonal parts of Eq. H2t)|) . We assume that the solution p{{p) can be expressed as p{{p) = 7to(p) +iZi{p) 
with 7^0 (p) and 7^i(p) as the lowest- and first-order terms of p{{p) in the A expansion. In the lowest-order of A, the 
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diagonal elements of the distribution function, TZofif^ij)), are determined by 

eE • Vppo^^ (p) = -i; [TZo] - [7^o] , (30) 

where the diagonal terms of Z^''^[7^o], depend only on the diagonal elements of 7^o(p), '^Omm(p)' ^^"^ given 
by Eqs. ljA.l|) and (|A.3|I . Since X^[7^o] is a higher-order term in the electron-impurity scattering, we can assume 
that X^[7^o] is much smaller than T^[TZo]. Hence, Eq. (|30|l can be solved as follows: we first solve Eq. (|30|l by ignoring 
[TZq] and then substitute the obtained solution into 2^ to calculate a correction caused by . Thus, we find that the 

solution of Eq. consists of two terms: 7?.o^^(p) = "^omm^p) ^O/i^'-P-'' ^^^^^ ^o^im'-P) ^^'^ ^0/i/i(p)' respectively, 
determined by 

eE.VpPo^^(p) = -i;[7^g], (31) 

and 

J^[7^S]-KX^[7^§]=0. (32) 

The solution of Eq. (jHlfl . 7?.g^^(p), depends on the momentum angle through sin(/)p and takes the form, 7^q^^(p) — 
eE^Q^{p) sin (pp, with the functions ^q^{p) given by 

t ^ _ OnF(£^(p)) [^lpp^^2fit,) dp +^3mm dp,, , s 

where is given by equation e^(pp) = £p(p) and the various relaxation times r^^iy (i = 1..3, /i, = 1, 2) are defined 
as 

— =2™,^|l/(p-k)pA,,,,(p,k), (34) 

with Ai^^(p,k) = i[l - cos((?!)p - 0k)]a+(p,k)A^^, A2^^.(p,k) = ia_(p,k)A^j,, A3^^(p - k) = i cos(0p - 
0k)a_(p,k)Ap^ and a±(p,k) = [XpXk ± ± a^fcpcos((/)p — (j)k))/XpXk- Substituting the result for 7?.g^^(p) 
into Z^, we find that the solution of Eg. I|32|l . Ttg^^ip), relates to a cosine function of the momentum angle: 
''^Omm(p) = e-B$g^(p) cos(/)p. %^ip) is given by 

/X _ ('^l/p + '^2jl,i) ^Op.iP) + Ts^fl^OjliPp.) .„p.x 
"^OplW - ~ / -1 , -I \ (-1 , ^-M _ ^-l ^-1 ' 



with 

'CS^(P) = TT^a^Mn, ^Vp_kVk_qyq-p-^Sin(,^k 
k,q ^P'^*^ 

x{c^+(p,q,k)A^^$g^,(p) +c^_(p,q,k)A^p$g^,(p^)}, (36) 

and c^±(p, q, k) defined by Eq. (|A.5|I . 

The off-diagonal 7^o(p) elements, 7?.o;ip(p), are obtained by substituting 7?.g^^(p) into the off-diagonal parts of 
Eq. 

*[iIo,7to(p)] = -5"[7t^]. (37) 

Here, the effect of 7^o/j/j(p) ^he off-diagonal elements was ignored. We find that 7^omm(p) depends on the momentum 
angle not only through sin(/)p but also through cos(/)p. However, as discussed above, we are interested only in 
the part of Re[(7?.o)i2(p)], which depends on 0p through cos0p, and the part of Im[(7?,o)i2(p)], involving sin^p: 
Re[(7?,o)i2(p)] = Co(p)cos0p + ... and Im[(7^o)i2(p)] = Co (p) sin 0p + with ^^(p) and Co(p) taking the forms, 

Co'(p) = E inP-k)pIm[g,(p,k)](-l)^{A^^ci>g^(p)[l-cos(0p-0k)] 

-A^p[ci>g^(p) - $g_(p^) cos(0p - 0k)]}, (38) 
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and 



4A^ 

k/i=l,2 

X {A^^$g^(p)sin(0p - <?!>k) - A^^$g^(p^)sin(0p - 0k)} , (39) 

and ga{p, k) defined as ga{p, k) = [afcAp sin((?!)p - (/)k) + iaM{p - fc cos((/)p - 0k))]/(A/cAp). 

To obtain tiie first-order term of the distribution function in tfie A expansion, 7?,i(p), we substitute 7?.g^^(p) into 
the diagonal components of the scattering term, and X^, as well as its off-diagonal component S''. We find 
that X^[7^q] depends on the angle of momentum through a sine function, while X^[7^g] relates to a cosine function: 
Ki'^o] = ^it^ip) sin0p and Xf^[7^g] = Cl^{p) cos^p with Cl'^ given by 

Cl^ip) - 2^n,Aa2^|y(p-k)|2:^sin2(0k-</>p) 

x{[l - cos(0k - 0p)]$L(p)A^^ - [$g„(p) - n-Jp^) cos(0k - '/'p)]AM/i} (40) 



and 



^l^^ip) = Tr^AAfn,^ Vp_kVi-qyq_pY^sin(0k 



k,q ^f^'^ 

X {d^+(p, q, k)A^^$g^(p) + (p, q, k) A^;,$g^(p^)}. (41) 

Here, o?^±(p,q, k) are defined by Eq. (|A.6|) . From the diagonal parts of Eq. I|26f) in the first order of A, Z°[7?.i] + 
2■^[7^g] +I^[no] = 0, it follows that 7^l^^(p) can be written as 

7^l^^(p) ^ e^$^^(p)sin0p -|-eS$J^(p)cos0p, (42) 

with $i^(p) determined by 

^s,cf , _ {'^IjifL + '^2p1») ^1^. (P) + '^Si^ti^lp (Pm) 

/ -1 , -1 W -1 , -M _ -1 -1 ■ ^'i^^ 



The off-diagonal elements of 7^i(p) are obtained from the off-diagonal parts of Eq. H26|) in the first order of A: 
j[J?o, '^i(p)] = ^'5^[72-o]- According to the definitions, Re[(7^i)i2(p)] — ^{{p) cos<j)p + ■■■ and Im[(7^i)i2(p)] = 
{p) sint^p + £,{{p) and Ci (p) can be written as 

Cf(p) = E inP-k)pIm[.g,(p,k)](-l)'^{A,,^<I>g^(p)[l-cos(0p-0k)] 

P kM=l,2 

-A^p[cl>g^(p) - <i>g^(p^) cos((/)p - 0k)]}, (44) 

and 

e((p) = ^fP E inP-k)pRe[g,(p,k)](-l)^sin(0p-0k){A^^<i>g^(p)-A,,<i>g^(p^)}. (45) 



kM=l,2 



The other component of the kinetic equation, Eq. H27|) . can be solved easily. The solution (p) has null diagonal 

'p + iC"{p) sin0p 



elements. Its off-diagonal elements can be written as (pi^)i2(p) — ^^^{p) cos (j)p + \p) sin^p with ^^^(p) and C^^(p) 



defined as 



and 



e'ip) = ^ (1 - [n^eiip)] - nAe2{p)]} (46) 



C"(P) = -^^{nF[siip)] - nF[e2ip)]}. (47) 
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C. Anomalous Hall current 

We first analyze the component of the anomalous Hall current, J^, that is associated with the current operator 
term arising from the free-electron Hamiltonian and is a sum of contributions from the diagonal and off-diagonal 
elements of the distribution function: = + ,P° . We know that J^^ is a component of the skew-scattering 
AHC. Considering Eq. (|13|l . it is obvious that the nonvanishing J^*^ comes from the diagonal terms of the distribution 
function, 7^op^(p) and 7^^^^(p), which depend on momentum angle through the cosine function. Thus, can be 
written as J^J^ = J^~^ + J^~^ , where J^~^ and J^^~^ , respectively, are associated with the lowest- and first-order 
terms of the distribution function, '^oaim(p) ^^'^ ^ip^(p)- 

^e'EY,{(^ + (-1)''^) Pcos' ^Pn.ip)] , (48) 

p,/j, ^ ^ ^ ^ 

= e^iS^ j ("i- + (-1)^^) pcos^ ^p^Im] ■ (49) 

Note that both J^^~^ and J^^~^ are proportional to the inverse of the impurity-density, i.e. (n^)"^, appearing when 
electron-impurity scattering is considered up to the second-Born approximation. 

Since both the distribution terms /5^(p) and /5^^(p) have nonvanishing off-diagonal elements, the contribution to 

the anomalous Hall current from off-diagonal elements of p(p), J^°, can be expressed as — J^°| + J^°| , where 
J^°|^ and J^\^^ arise from /5(2(p) and /5(|(p), respectively, and take the forms 

Ji°\' = 2e'Ej2 [Co(p) + cos' ^ + a [Co'(p) + (l{p)] sin' i 

and 



(50) 



= ^"^^ E I? (1 - V) {uFlsiip)] nF[£2(p)]}. (51) 

^ p \ 

From Eq. H5UI) we see that J^°|^ is independent of the impurity density, due to the n^-independence of and Ci- 

f 1 1 

However, it is due to disorder and relates to longitudinal transport. It is obvious that involves the derivative of 

f I / . 

the equilibrium distribution function, i.e. dn-p{uj) / duj. This implies that J °\ arises only from electron states in the 

vicinity of the Fermi surface, or, in other words, from electron states involved in longitudinal transport. Physically, 

the electrons participating in transport experience impurity scattering, producing diagonal p((p) elements of order of 

n,~^. Moreover, the scattering of these perturbed electrons by impurities also gives rise to an interband polarization, 

which eliminates dependence on the impurity density within the diffusive regime. 

fill. 

However, the anomalous Hall current J °\ is a function of the entire unperturbed equilibrium distribution, n-piuj), 

not just of its derivative, dnF{uj)/duj, at the Fermi surface. This indicates that has contributions from all 

f ill . . 

electron states below the Fermi sea. Obviously, is independent of any electron- impurity scattering and relates 

to the driving terms, one of which is just the interband electric dipole moment, while the other one arises from the 
SO coupling directly induced by the driving electric field. 

From Eq. (|10|l . it is obvious that to determine the first-order term of J™^ in the A-power expansion, one has to 
deal with the function < V'^p^/'^k > in the lowest order of A. We find that this lowest-order term of < ip^pipfik > can 
be evaluated from the kinetic equation for the distribution function p{p). To show this, we start with a Heisenberg 
equation for the operator V'iJp^A'P- 

^r^^V^^p^MP = [H, ^tj,p] = -»eE • Vp[^+p^^p] + [Ho + He, V^+pV^^p] + I%{p), (52) 
where, /^^(p) = [/fimp, ^itp^Mpl- lowest order of A, /^^(p) takes the form, 

^;.(P) « E ^P-k[e'^'-(P-'^)^+k4p - e'^'-^'^-P^V'i^Mk]. (53) 
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Multiplying both sides of Eq. H52() by eimnPmO'^^ and taking the summation over /i and we get 

J2 ^lmn[Pma:^ < /^.(p) >] « ^ e'""' ' ^""^^ ^k-p£i™ (fc™ - Pm)[< i'+J^W > J , (54) 
Pi/^^ p,k,/ii' 

with < ... > denoting a statistical average. Obviously, the right-hand side of Eq. (|54|) is just the AHC component J™''. 
On the other hand, taking the statistical average of Eq. (|52|l reduces it to the kinetic equation for p^j, (p) =< tp^ptp^p >■ 
Hence, from Eqs. (|52|l and 154(1 . it follows that J™^ can be written in the spin basis as 



J™P = -Ae^£,™„p„Tr{a" (eE • Vpp(p) + i[7?o, /5(p)]) } , (55) 
p 



where the contribution associated with He is ignored because it is of higher order in A. J™^ can also be determined 
in the helicity basis by means of 

J™P = -Ae ^irnnPrr^Tr [u+c7^Up [cE • Vp/5(p) - eE • [p(p), U+VpUp] + t[Ho, p(p)]] } . (56) 
p 

In the linear response regime, Eq. 1)56(1 reduces to 

JT" - ~XeJ2^{eEMV,^[{poUip)-{po)22ip)] 
p ^ 

_eEa^^^ sin</)p[(po)ii(p) - (/5o)22(p)] + 4apAp[Co (p) + C"(p)] sin(^p 

= -Ae^] {nF[ei{p)] - nj. [e2(p)]} + 4app,[Co'(p) + C" (p)] sin^pj , (57) 

To obtain the last equality in Eq. 1(571) . the momentum integral with integrand involving the derivative with respect 
to Py was performed by parts. Note that Sinitsyn, et al. recently investigated this component of the anomalous Hall 
current in the absence of Rashba SO coupling by analyzing the effect of impurity scattering on the coordinate shift 
of the electron wave-packet^ 

arises from SO coupling directly induced by the driving electric field. Considered to linear order in the electric 
field, it takes the form 

Jf = Xe'EY^iuFleiip)] - nj. [e2(p)]} . (58) 
p 

Obviously, this contribution to the anomalous Hall current is independent of any electron-impurity scattering. If only 
one of the parameters-the Rashba SO coupling constant or the magnetization-is zero and other remains finite, 
doesn't vanish. In contrast to this, Jt^^ and as well as J™p reduce to zero for just one of them vanishing, null a or 
M. Note that in all previous studies, the contribution to anomalous Hall current from the SO term due to the driving 
electric field has been ignored. In the following numerical calculation, we will show that plays an important role, 
especially in the low magnetization regime. 

Thus, after all components of AHC are determined, the total anomalous Hall current can be obtained from Eq. 
We define the total anomalous Hall conductivity as a^y — Jx/E. Obviously, a^y can be written as 

With a% = .rJlE = a% + ai-P -t- af^ and a- = Jf/E = af-^ + a^-^ . Here, the quantities, a%, a^, a^^S 
and are defined as a^P-^ = J^P'^/i?, < = j'^jE, af-^ = Jf-^IE and a--^ = J^r^ IE. 

It should be noted that, in our study, the diagonal part of the distribution function, 7?.q(p), which is involved in 
all disorder- related components of anomalous Hall current, was evaluated in the self-consistent Born approximation. 
This implies that our results correspond to that obtained in the Kubo formalism by considering the "ladder-sum" 
vertex corrections to the bubble diagrams. 
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III. RESULTS AND DISCUSSIONS 



We have carried out a numerical calculation to investigate the anomalous Hall effect in a InSb/AlInSb quantum 
well with Rashba SO coupling. Such a system was recently examined experimentally^ It is well known that the 
InSb semiconductor is a good material for AHE observation because its band gap, Eq = 0.235 eV, spin-orbit splitting, 
^SO = 0.81 eV, and P — 9.63 eV-A result in a pronounced spin-orbit coupling constant A = 5.31 nm^ (for GaAs, 
A = 0.053nm^)iSi Also, the large g-factor, g = —51.4, may lead to a remarkably large magnetization. In our 
calculation, the static dielectric constant, k, and the effective mass of InSb, m* , are chosen to be k = 17.54 and 
m* = 0.0135mo with mo as the free electron mass. The width of the InSb/AlInSb quantum well is assumed to be 
a = 20 nm and the density of electrons is taken as iVg = 1 x 10^^ m~^. We consider an attractive interaction between 
the electrons and the background impurities in the quantum wells (the attractive and repulsive interactions lead to 
differing anomalous Hall effects because their contributions to AHC in the second Born approximation have opposite 
signs). Note that we have also estimated the effect of scattering of electrons by remote impurities on AHE, finding 
that it is relatively small and can be ignored. Thus, the scattering potential Vq can be written as2i 



Vq^U{q)F{q)/K{q,0), 
with U{q) — —e'^/{2eoKq) and the form factor F{q) determined by {u = qa) 



(60) 



(47r2 -I- u'^)u 



u 1 — cxp(— m) 
47r^ u 



(61) 



K(q, 0) is a static dielectric function in random phase approximation and can be written as 



with qs = 7Ti*e^/(27reoK) and H{q) given by2^ 
^ 1 — exp(— w) u 



H{q) 



K{q,0) = l + ^Hiq), 



1 — exp(— u) , o o 2 
^ -(u^-47r^) + - 



(ti2 + 47r2)2 



1 - 



1 — cxp(— u) 



(62) 



(63) 



Here, the effect of the Rashba SO coupling on the screening of is ignored. Further, to determine the impurity 
density, we assume that for AI — and a = the electron-impurity scattering results in an electron mobility 
/xo = 5m2/Vs. 



A. Anomalous Hall effect in a InSb/AlInSb quantum well without Rashba SO coupling 

We first analyzed the anomalous Hall effect in the absence of Rashba SO interaction. In this case, the lowest order 
component of in A, J^^^^, vanishes, while the component J^^~^ is nonvanishing and reduces to 

= e'Ej2£^ cos' c^p^.ip), (64) 

P,A« 

with ^i^{p) determined by 

^uAp) = (-l)^47r2n,ri^^A$g^,(p) ^ Vp-kV"k-qVq_ppfc sin((?!)k - 0p)(5(e^p - e^fe)(5(e^p - e^,) 

x[pksm{(f>i^ - (j)p) + qpsm{(l)p ~ 0q) - qksm{(f>i^ - 0q)], (65) 

and ^Qfj^ip) — — Ti^p ^"a}'^(p) ^^ ^^dp''^ • ^i^*^^ the contributions to AHC from the off-diagonal elements of the distribution 

function, ,P° = J^°\^ J^x\^^ ^ vanish, the side-jump AHC involves only the components J™^ and J^: — J^^^ + J^- 
Here, J""^ and are equal to each other and take the form 

JT" = = Xe'Ej2{nF[ei{p)] - np[e2{p)]} . (66) 
p 

In Fig. 2, we plot the calculated total anomalous Hall conductivity axy = cr^y + cr^^y, and its component trj'^ = 
(T™P + '^xv functions of magnetization, M. With increasing magnetization, axy and cr|' increase linearly. A 
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M (meV) 



FIG. 2: Magnetization dependencies of Oxy and a^y in a InSb/AlInSb quantum well without Rashba spin-orbit coupling. The 
width of the quantum well a is a = 20 nm. The electron density is TVe = 1 x 10^^ m~^. The lattice temperature is T = OK and 
the mobility in the absence of magnetization is — 5m^/Vs. 

comparison between Uxy and (j^]y indicates that, for the given /xq = 5m^/Vs, both the contributions from side jump 
and skew scattering are of the same order of magnitude. Note that, notwithstanding the large spin-orbit coupling 
constant A, the anomalous Hall conductivity is still much smaller than the ordinary one: the ordinary Hall conductivity 
is 34.6e^//i for a magnetic field B — 0.34T (in the InSb/AlInSb quantum well with g ~ —51.4, this magnetic field 
corresponds to a magnetization M = 1 mcV). 

From Eqs. (|65|l and (|66|l we see that in the absence of both the Rashba SO coupling and magnetization, the 
contributions to anomalous Hall current from electrons with opposite spins (or hclicities) have opposite signs. As 
a result, the total anomalous charge Hall current vanishes. However, there is a nonvanishing spin Hall current 
since electrons with opposite spins move toward opposite sides of the sample. We estimate the spin-Hall current 
in the studied InSb/AlInSb quantum well for a = and M = 0, finding that the spin-Hall mobility /i^if defined 
in Ref.im is ^sH — 0.013 m^/Vs. (In contrast to this, in a GaAs/AlGaAs quantum well with charge mobility 
/io = 0.6m^/Vs, the total spin-Hall mobility is jisH — —2.0 x 10~^m^/Vs, and the contributions from side-jump 
and skew scattering, respectively, are /ig'^ = —1.6 x lO^^m^/Vs and /i^lf — 1.4 x 10~'*m^/Vs. They are of the same 
order of magnitude as the spin-Hall mobilities in bulk n-doped GaAs: in a bulk GaAs with the same /ioi Engel et 
al^ found fil^ ^ -1.6 X IQ-^mVVs and ^'/^ = 3.5 x lO-^m^/Vs.) 

B. Anomalous Hall effect in a Rashba InSb/AlInSb quantum well 

We have also calculated the anomalous Hall conductivity in a InSb/AlInSb quantum well with Rashba spin-orbit 
interaction. In the case of nonvanishing a, one has to consider not only (J^y~^ and o-"y^'^, but also the anomalous 

Hall conductivities, (t^° = '^xyl^ + ""Syl^^; ^^'^ '^Ty'^- The results are plotted in Figs. 3 and 4. 

In Fig. 3, we plot the total anomalous Hall conductivity, a^y, and the skew scattering and side-jump Hall conduc- 
tivities, a^y and a^^y, as well as their components, ct^t^^, fJ^xyi '^xy: cr™P, as functions of magnetization M for 
various Rashba spin-orbit coupling constants. We find that, as magnetization increases, the total anomalous Hall 
conductivity axy increases for a = 5, 10, and 20meV-nm, but it decreases for a — 50meV-nm. 

Such complicated behavior of the magnetization dependence of axy arises from competition of the side-jump and 
skew-scattering contributions to anomalous Hall conductivity. From Figs. 3(b) and 3(c) it is obvious that a^y varies 
monotonically with the magnetization, but in the magnetization dependence of a^'y there is always a small dip (when 
a = 50meV-nm, this dip is shifted out of the studied magnetization range). In a^y, cr|y~^ is dominant for small 
a, leading to an increase of a^y with increasing M. However, when a > 30meV-nm, cr^y'^ is important: a^y 
becomes negative and its magnetization dependence exhibits a decrease with increasing magnetization. Among the 
contributions to a^^y, a^°y is dominant for large magnetization, while a^y is important for small M and results in a 
nonvanishing a^^y when a ^ Q but M — Q. 

In Fig. 4, we plot Uxy, f^y, and cr^Jy as functions of the spin-orbit coupling constant. We find that, as a increases, 
Uxy increases monotonically for M = OmeV, while for M = 0.5 — 2meV, cTxy first decreases and then increases. It is 
also evident from Figs. 4(b) and 4(c) that, with increasing Rashba spin-orbit coupling constant a, increases while 
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FIG. 3: Magnetization dependencies of (a) the total anomalous Hall conductivity, (Jxy, (b) the skew-scattering and (c) side-jump 
Hall conductivities, as well as their components, (d) cr^xy^^ ifi) (f) '^xy i and (g) al^^ in a InSb/AlInSb quantum well with 
different Rashba SO couplings: a = 5, 10, 20, and 50meV-nm. The other parameters are the same as in Fig. 2. 



(T^y decreases. Since the rate of increase or decrease depends on M, the a dependence of the total anomalous Hall 
conductivity behaves differently for different magnetizations. 



IV. CONCLUSIONS 



We have employed a kinetic equation approach to investigate the anomalous Hall effect in Rashba 2D electron 
systems based on narrow band semiconductors. The Rashba SO coupling was considered nonperturbatively, while the 
extrinsic spin-orbit interaction and the SO coupling directly induced by an external driving electric field were taken 
into account in the first order of the coupling constant. Considering electron-impurity scattering up to the second- 
Born approximation, we found that the various components of the anomalous Hall current can fit into two classes: the 
side-jump and skew scattering anomalous Hall currents. The side-jump anomalous Hall current involves contributions 
not only from the extrinsic SO coupling, but also from SO coupling directly induced by the driving electric field. 
It also contains a component which arises from Rashba SO coupling and relates to the off-diagonal elements of the 
hclicity-basis distribution function. The skew-scattering AHE arises from the anisotropy of the diagonal elements 
of the distribution function, and it is a result of the Rashba and extrinsic SO interactions. We also performed a 
numerical calculation to investigate the anomalous Hall effect in a InSb/AlInSb quantum well. We found that the 
contributions to anomalous Hall conductivity from both the side-jump and skew scattering terms are of the same 
order of magnitude, leading to complicated dependencies of the total anomalous Hall conductivity on magnetization 
and on the Rashba spin-orbit coupling constant. It is also clear that the component arising from the SO coupling due 
to the driving electric field is dominant for small magnetization. 
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FIG. 4: The dependencies of (a) the total anomalous Hall conductivity CTxy, and (b) the skew-scattering and (c) the side-jump 
anomalous Hall conductivities, a^y and a^'y, on the Rashba spin-orbit coupling constant for different magnetizations: M = 0, 
0.5, 1.0, 1.5, and 2.0 meV. The other parameters are the same as in Fig. 2. 
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APPENDIX: EXPLICIT FORM OF /^^^ 



The diagonal elements of can be written as sums of four terms: = I^[pi] +'^^[pi] +^^[pi] + 

I^IPi] {[pi] denotes that J^.^.^.d depend on the specific form of the distribution function p{). I^[pi] and I^[pi] are 
the terms of the first Born approximation and can be written as 



k 

+a_(p,k)[(p()^^(p) - (p()pp(k)]A^p} , 



(A.l) 
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and 



Il[p{] = 27rn,Aa2^|np-k)| 



XpXk 



(A.2) 



with A^y = (5[ep(p) — eiy(k)] and a±(p, k) = (ApAfc ± ± a^fcpcos(0p — <?!'k))/ApAfc. I^[/5(] and I^[/3f] are the terms 
of the second Born approximation. I^[pi] is exphcitly independent of the spin-orbit coupling constant A and takes 



the form 



k,q ^P^*= 

X {c/.+ (p,q,k)A^^(p()^^(k) +c^_(p,q,k)A^p(p()^p(k)} 



while I^[pi] is linear in A and can be written as 

k,q 



ApAfe 



X {d^+(p,q,k)A^^(p{)^^(k) +d^_(p,q,k)A^^(p{)^;i(k)} . 
The parameters c^±(p,k, q) and d^±(p, k, q) are defined as 



c,.±(p,q, k) 



A„ 



-C(p,q,k)[5(e^(p) -£^(q)) - 5(e^(p) -£^(q))], 



(A.3) 



(A.4) 



(A.5) 



with 



and 



ci^±(p,k,q) = -{I?™±(p,q,k)[5(£^(p)-£^(q))-5(£^(p)-£^(q))] 
+P^±(p,q,k)[5(£^(p) - s^id)) + 6{s^{p) - e^(q))]} , 



C(p,q, k) = pqsm{4'p - 0q) - qksin{4>k - 4>q) +pfcsin(0k - M, 

X>™±(p,q,k) = {-1Y+\±M^ + XkXp)\pksm{(j)i, - (t)p) + qpsm{(f>p - ^q) - qksm{(l)i, - ^q)] 

-|-(_1)/'+1q,2 {^pksm{(f>i^ - (j)p)\pqcos{(j)p - (^iq) + qkcos{(j)ci - ^k) - kpcos{(j)u - (pp)] 
+gpsin(0p - (f)q)[-pq cos{(j)p - (j)^) + qkcos{(j)q - ^k) + fcpcos((/>k - (pp)] 
+qksm{(j)i^ - (/)q)[-pgcos((/)p - ^q) + qkcos{(j)q - t^ik) - kp cos{(pu - <l)p)]} , 



P^±(p,q,k) = (-l)''+^Ag(Afc ± Ap)[pfcsin((?ik - ^p) + qpsm{(l)p - 0q) -qksm{(f)k - 0q)]. 



(A.6) 
(A.7) 



(A.8) 



(A.9) 



Since the off-diagonal elements of the collision term /^^^ are simply related by /^2^ = —[-^12^' suffices to consider 
the element I^^' . In the first-Born approximation, 7^2^ can be expressed as a sum of two terms: ^2^ = S''[pi]+S''[pi] 
with <S°[/5(] and »S''[/3(] as the terms in the lowest- and first-order of A, respectively, and determined by 

^"'"[Pi] = yEl^(P-k)l'5a,6(P>k){[(GS.)ll(^2p)-(GS,)22(£2p)] (p()22(p) 
k 

- [(GSp)22(£2fc) - (GSp)ll(e2fe)] (p()22(k) - [(GSfe)ll(£lp) - (Go%)22(£lp)] (/5{)ii(p) 

+ [(Ggp)22(£ife) - (GSp)ii(£u)] (/5()ii(k)} . (A.10) 
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Here, ga(p,k) = [afcAp sin(<^p - (pu) + iaM{p - fccos(0p - 0k))]/(AfeAp) and gb{p,'k) = -2sin(0p - (j)k.)[iM sin{(j)p 
M + ApCos((/)p - (j)i,j)]XaPp/{XkXp)- 
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